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We revisit the topic of atomic center of mass motion of a 
three level atom Raman coupled strongly to an external laser 
field and the quantum field of a high Q optical cavity. We 
focus on the motion related nonadiabatic effects of the atomic 
internal dynamics and provide a quantitative answer to the 
validity regime for the application of the motional insensitive 
dark state as recently suggested in Ref. [Phys. Rev. A 67, 
032305 (2003)]. 
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I. INTRODUCTION 

The development of quantum information science and 
technology carries the potential of revolutionary impact 
on many aspects of our society, as evidenced already 
by the applications in quantum cryptography, quantum 
communication, and rudimentary quantum computing. 
Among the physical systems being investigated, high Q 
optical cavities coupled with trapped atomic qubits rep- 
resent a paradigm for this burgeoning field. In addition 
to their demonstrated abilities for controlled (and co- 
herent) quantum dynamics of both atomic and/or cav- 
ity photonic qubits, cavity QED systems are unique be- 
cause they represent a proto-type enabling technology 
for the coherent inter-conversion of quantum information 
encoded in material qubits or flying photonic qubits for 
propagation to far away places. 

Despite much effort and spectacular advances from 
several groups in recent years [1-3], high fidelity deter- 
ministic logic operations even at the level of two qubits 
remain elusive in cavity QED based systems. Among the 
factors as commonly attributed to being significant road 
blocks, the localization of atomic motional wave packet 
is perhaps the most demanding. In nearly all quantum 
computing protocols of atoms coupled to a high Q cavity 
field, it is essential to reach the so-called strong coupling 
limit, where the coherent coupling of an atom with the 
near resonant cavity mode g has to be much larger than 
both the cavity decay rate k (one side) and the atomic 
spontaneous emission rate 7, i.e. g ^> k and 3^7. 
Since g is inversely proportional to the mode volume 
of the cavity, it typically points to small cavities in the 
Fabry-Perot arrangement, where the cavity mode is that 
of a standing wave given by 

9(r)= 9oX(r), 



with p = \/ x 1 + y 2 the transverse (polar) coordinate 
measured from the cylindrically symmetric cavity axis 
along z direction. The typical geometries have the mode 
waist w(z) = Wq^/I + z 2 /z 2 much larger than the cavity 
wavelength A, where z — ttWq/X is the Rayleigh range. 
Unless each individual atomic motion is localized to much 
less than the resonant wave length A, i.e. in the so-called 
Lamb-Dicke limit (LDL), this position dependent uncer- 
tainty of coupling strength g{r) generally spoils the quan- 
tum coherence, and prevents high fidelity quantum logic 
operations. This challenging limit is not so- far under 
complete experimental control, it is especially problem- 
atic for optical cavity QED systems. 




FIG. 1. Illustration of the Raman coupling. 

Recently, two independent groups [4,5] have noticed 
an interesting scenario where the above undesirable po- 
sition dependence of g(r) can be largely overcome with 
the use of a so-called dark state, when the classical Ra- 
man laser field is assumed to have the same spatial de- 
pendence as the quantum field g{r). For a large class 
of quantum computing protocols based on atomic cavity 
QED, the building block consists of a three level A- type 
atom of stable ground states |po) and \g\) that couple 
to an excited state |e). Such an arrangement allows for 
a coherent mapping of an atomic qubit (a\go) + /3\gi}) 
into the photonic coherence of the cavity. In the most 
publicized version as originally suggested [6] , it is as- 
sumed that a classical laser field and the quantum cavity 
field establishes a two-photon matched Raman resonance 
between \g±) <-> |e) and |<7o) By denoting the 

Rabi-frequency of the classical laser field as tt, the dipole 
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coupling in the interaction picture can be summarized as 

H = -hA\e)(e\ 

+hn\e)( gi \ + hg\e)(g a \c + h.c., (2) 

where c is the annihilation operator for the near resonant 
cavity photon mode used for Raman coupling, and the 
common detuning A = lol — w egi = u) — Lo ega . It is easy 
to check that the following "dark state" 

l £) > = Trwrm (glgi.°)-%o,i)), (3) 
vwr + l^l 2 

is an eigenstate of the Hamiltonian Eq. (2) with a zero 
eigenvalue, where |0)c and |l)c denote the Fock state of 
or 1 cavity photons respectively. It is dark, or, immune 
to atomic spontaneous emission because it contains no 
atomic excited state. By engineering a counter-intuitive 
pulse sequence as in the STIRAP (stimulated Raman adi- 
abatic passage), and assume the atom + cavity system 
to adiabatically follow the above dark state, it leads to a 
highly efficient protocol for converting the atomic qubit 
state into a photonic superposition according to [6-8] 

(a\g ) + P\gi)) ® |0)c -> \gi) ® (a|0) + 0\l)) c . (4) 

When atomic motion is considered, the position de- 
pendence of g{r) generally leads to a loss of coherence 
due to the potential entanglement between the motion 
and the atomic internal state as well as the cavity pho- 
ton state. The idea of the motional insensitive protocol 
[4,5], assumes a classical laser field that has the same de- 
pendence as g(r). Following the notation of Duan et. at 
[5], this assumption amounts to 

0(f, t) = n (t)x(0 = r goa(t) X (r). (5) 

A simple arrangement involves choosing the pump and 
the cavity transitions \gi) <-> |e) and \g ) <-> \e) to cor- 
respond the left and right circular polarized component 
of the same cavity mode. A more flexible setup would 
involve the use of a different cavity mode such that near 
the cavity center, x{r) f° r the two modes remain almost 
matched as the two modes differ very little in their re- 
spective wavelengths [9]. 

The aim of this paper is to study in detail nonadia- 
batic effects of the above motion insensitive protocol. As 
was also noted in the Ref. [5] , it clearly becomes difficult 
to maintain adiabaticity when the atomic Raman cou- 
pling is too weak to affect the transfer, particularly near 
regions of small g(r) values. Furthermore, an atom re- 
maining in the dark state essentially experiences no light 
force from the combined fields of both the cavity mode 
and the external laser. This arguably leads to an upper 
limit on the atomic kinetic energy; the duration T for the 
STIRAP is determined by atomic internal state dynam- 
ics. Thus for an atom with a velocity of v a , during the 
STIRAP, it will move a distance v a T if it is to remain in 



the dark state. A larger v a leads simply to a large trav- 
elling distance. Nonadiabatic effects will arise if atoms 
were to travel far enough to cross nodal planes of \{f) 
(as we will see later this contradicts the discussion of a 
better satisfied adiabatic condition near nodal points as 
in Ref. [5]). The ideal operation of the motional insen- 
sitive protocol would require the use of trapped atoms, 
i.e., with atoms confined near regions of maximal g(f) 
[5] by an external force independent of the cavity or Ra- 
man field. The trap has to be strong enough to limit the 
atomic motion due to an amplitude smaller than half the 
cavity wavelength A/2, such that nodal crossing can be 
completely avoided. 

Another motivation for this study is the desire to un- 
derstand nodal crossing dynamics in general for the Ra- 
man configuration when the trapping provided by an ex- 
tra higher order cavity mode is absent as in earlier ex- 
periments [2] , where either the resonant cavity field or an 
external laser field provided confinement of atoms (not in 
the dark state). Based on the dressed energy levels of an 
atom coupled to both fields, we find that the dynamics 
of an atom crossing a nodal plane can be effectively de- 
scribed in terms of the celebrated Landau-Zener theory 
[10]. Yet, a somewhat puzzling situation arises according 
to Landau-Zener theory which prefers to have a larger ve- 
locity v a during the crossing in order to maintain in the 
initial (dark) state. Consistent with the paper of Duan 
et. al. [5], the motional state insensitive protocol works 
only in the limit when the atoms are trapped by yet an 
additional mechanism such that its motion is limited to 
a variation of g(r) within approximately a factor of 2. 
On the other hand, a large velocity tends to cause large 
amplitude motions, thus against the localization of the 
Lamb Dicke limit. We thus find it interesting to study 
the relevant Landau-Zener transitions in order to shed 
light on the motional effects of atoms in cavity QED. 

This paper is organized as follows. In section II, 
we formulate the model of our study and illustrate pa- 
rameter regimes of interests to current experimental ef- 
forts. Sections III and IV are devoted respectively to 
the study of the nonadiabatic level crossing in terms of 
a Landau-Zener transition dynamics and the comparison 
between numerical simulations and the approximate an- 
alytic Landau-Zener state transition formulae. We have 
developed an interesting analytic mapping (in the ab- 
sence of an external trap) of the atomic motion through a 
nodal point into a Landau-Zener level crossing dynamics. 
Within each of the above sections, we will study various 
limiting cases, mainly focusing on a simple model that 
involves a 1-dimensional motion along the cavity axis 
[11,12]. Finally we summarize and attempt to make some 
general conclusions in Sect. V. The appendices contain 
several technical points that may be useful for related 
studies. 
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II. FORMULATION 

In the descriptions to follow, we will assume both fields 
to be on resonance and take the atomic detuning A = 0. 
For the more general situation as shown in the appendix 
A with a nonzero but constant A ^ (position inde- 
pendent), we find it simply leads to formally identical 
results as discussed here for A = (see appendix A). 
When necessary, an external trap, assumed to be internal 
state independent is assumed to be available and centered 
around the maximum of g(f) to confine atomic motion 
to within the order of half the standing wave wavelength 
as in Ref. [5]. 

The Hamiltonian for atomic internal degrees of free- 
dom (on resonance A = 0) can be written as a 3 <g> 3 
matrix [5] in the basis {|e,0), \go, 1), |<7i,0)}, 



H 




(6) 



Within the semiclassical approximation for the atomic 
motion, the center of mass motion is described by p 2 /2M 
with M the atomic mass and p the atomic momentum. 
While for slow atoms, we may wish to include the as- 
sociated forces due to the atomic interaction with the 
spatially dependent laser (cavity) fields. We will first of 
all, consider the simple case of predetermined atomic mo- 
tion as corresponds to atoms staying in the dark state, 
i.e. we simply assume that atomic motion is not affected. 
We will also comparatively address the case of a trapped 
harmonic atom motion. We defer the inclusion of atomic 
dipole forces to a future investigation. With these as- 
sumptions, we find the three eigenvalues 



= 0, 

E ± = ± X [r(t)]g y/l + \r a{t)\*, 



(7) 



with E± depends 
Q = r a(t)g, \g\ 



on r. Here, 
= \xW)]9o\, and 



we have used 



r a(i)| 2 

X[r(t)]g > 0, and E+ < E_ when x[r(t)]g Q < 0. 



We note that E + > E_ when 




FIG. 2. Dressed energies of the coupled atom + cavity sys- 
tem at a particular radial location along the cavity axis z (at 
a fixed time t). When the atom moves through a nodal plane, 
nonadiabatic transfer of its internal state out of the dark state 
becomes a critical issue, which constitutes the main topic be- 
ing studied in this work. 

Figure 2 shows the atomic dressed state energies (7) 
along the cavity axis at a given radial location (A = 0) 
for a fixed time t when the external laser is tuned on. 
For an atom to cross a nodal plane, it has to be initially 
within a distance reachable within the duration of the 
external pump pulse a(t) with its initial velocity. 

In addition to dark state (3), the two other eigenstates 
are 

1 



B- , =-== (|fl)±|e,0», 



\B) 



oi<l> 



- : [\g o ,l)+r*a(t)\ gi ,0)}, (8) 



where e 1 ^ = g*/\g\ and will assumed to be unity later. 

Before writing down the Schrodinger equation, we ab- 
sorb the phase factors due to the adiabatic evolution 
along each of the above (time-dependent) eigenstate 



|*(t))= J2 C n (t)e- l Io E - Wdt ' /R \n), 



(9) 



n=0,+,— 



where we have used the shorthand notation of n = for 
the dark state \D) and n — ± for states \B±) respectively. 
The coefficients C„ are governed by 



C 



C+ = - 



e- l $o E+{t ' )dt ' /Tl {D\B + )C + 

e - l Io E -^ dt '/ h (D\B^C-, 

e +i Sl E+{t ' )dt,/n {B + \b)C 

e - l J>+^- E -^ dt '/ n (B_\B + )C, 



(10) 



We now assume that a(t) is a real parameter, which 
leads to 



(B±\B^ 



1 



(B\B) 



0, 



(D\B ± )= -j={D\B) 



v/2(l + |r a(t)| 2 )' 



(11) 



Denoting E± — ±e, we end up with the simplified form 
ofEqs. (10) 

C7 = -e-*Jo e{t " )dt ' /R K(t)C+ - e l Io e{t ' )dt ' /h K(t)C-, 
C + = e +t tt £(t ' )dt ' /h K*{t)C^ 

C-=e~ l fo e ^ dt '/ h K*(t)C , (12) 
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with 



K(t) 



r* a{t) 



V2(l + |r a(t)| 2 ) 
(D\B±). 



(13) 



These equations can be numerically integrated us- 
ing standard algorithms to investigate nonadiabatic level 
crossings. Before attempting an analytical understand- 
ing of the level crossing dynamics near the nodes of the 
cavity mode function x( 7? ) m the next section, we first 
consider here typical regimes of system parameters. 



A. parameters 

We will use Cs as a proto-type atom for the estimation 
of various atomic parameters. The resonant transition 
between the 6 2 Si/ 2 and 6 2 P 3 / 2 , (unclear spin I = 7/2), 
occurs at around A = 852.35 (nm), hu = 1.455 (eV) and 
excited state lifetime (I/7) = 30.70 (ns), or 7 = (2tt)5.18 
(MHz). The recoil frequency of resonant transition is 
about lor = hk 2 /(2M) — 2.07 (kHz), corresponding to a 
temperature of 0.198 (/uK), or a recoil velocity of 0.352 
(cm/s). The wave length of Afort = 936 (nm) used for 
the dipole trapping of atoms in Ref. [5] corresponds to a 
higher order cavity field. 

The vacuum Rabi coupling between the cavity field and 
the atom is go = (2*7r)50 (MHz) as in the recent CalTech 
experiment [1]. For order of magnitude estimates, the 
ratio r = fl/[ga(t)] can be taken as unity. 

According to Sect. III. (A) of Ref. [5], we take the 
adiabatic parameter a(t) as a Gaussian function 



a(t) = ao exp 



(*■ 



t 2 



(14) 



with the total protocol for state transfer being approxi- 
mated as T = 2T . The amplitude ao is assumed to be 30 
and the width tw is assumed to be Xq/3. For the efficient 
operation of the quantum state transfer protocol utiliz- 
ing the dark state adiabatic passage, the excited state 
atomic life time oc I/7 essentially sets the time scale, for 
this reason we take tw ~ I/7, or To <~ IO/7. 

Assuming that atoms are trapped inside a single well 
of the standing wave cavity field, or of the wells of the 
dressed energy E±(r), we can estimate the oscillation 
frequency inside according to 



9oX(P,z- A/4)=.g exp 
= .9o exp 
w g exp 



w 2 (z - 


-A/4) J 


P 2 1 


w 2 (z - 


-A/4) J 


P 2 ' 


(kz) 2 


w 2 (z)_ 


1 2 



sm(kz — 7r/2) 
[2sin 2 (^)-l] 



1], 



(15) 



near the axis center where the well is deepest, which gives 
the strongest axial oscillation and radial oscillation ac- 
cording to 



X -Mu? z z 2 = -Kg k 2 z 2 , 



Mw 2 p 2 = hg 2? v 



i.e. we obtain 



uj z = ^2go(hk 2 /2M), 



uj p = 2^g (hk 2 /2M) ■ 



1 



(16) 
(17) 

kw(z) ' (18) 

Given the additional enhancement due to Q(r) in E±, 
and use the estimated parameters as outlined above, we 
then take u z ~ (27r)500 (kHz), and assume a fundamen- 
tal cavity mode waist of w(0) = 30A, we end up with 
uo p <~ (27r)2.65 (kHz), of the same orders of magnitude 
as in Ref. [1] of the optical trap from a higher order cav- 
ity mode. Of course, these estimates are valid only for 
atomic motion near the bottom of the trap. For signif- 
icantly higher atomic energies, as for instance in the re- 
cent experiment [1], where atomic kinetic energy is of the 
order of (27r)20 (MHz), or about half of the actual poten- 
tial barriers at about (27r)50 (MHz), and about 40(hw z ) 
and 7500 (hu) p ), simple harmonic motion cannot be as- 
sumed. 

Within the regimes of these parameters, atomic center 
of mass motion is well approximated by classical dynam- 
ics in a conservative potential. We can also estimate the 
maximal velocity of these atoms when trapped in the 
above single well potential, 

±Mv 2 zM (Kk/M) 2 = A0(hu; z ) 

~ \Mv 2 pM {hk/Mf - 7500(7iwp), (19) 

where we have expressed velocities in atomic recoil units. 
Thus we find 



V Z M ~ V p M 



= 100(v R ) w 35 (cm/s). 



(20) 

This leads to the assumption of atomic velocities with 
the following choices for numerical simulations; v z w 1 
(m/s) <~ 0.036 (A7) (10 times more kinetic energy), v z w 
0.35 (m/s), and v z w 0.1 (m/s) (10 times less kinetic en- 
ergy). With To <~ IO/7, the respective distances a typical 
atom travels during the state transfer protocols become 
0.03 — 0.3 (/im), which is a significant fraction of A/4, 
half the distance between the nearest nodal planes. It 
is important to emphasize that in this limit which cor- 
responds to the recent experiment [1], atomic kinetic en- 
ergy is much higher than its single photon recoil energy, 
and the atom's motional quantum state is much higher 
than the ground state of each trapped well. This lends 
strong support to our assumption of using a constant 
(predetermined) atomic trajectory in studying the level 
crossing dynamics. When an additional cavity field is 
used to confine atoms, optical dipole force is the reason 
for atoms to turn around at classical turning points. 
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B. Qualitative picture of the failure of adiabaticity 

To maintain adiabaticity during atomic motion, the 
system must satisfy [13] 



(n{t)\d t \m{t)) 



E n (t)-E m (t) 



«1, 



(21) 



for all adiabatic energy levels. The above notation ap- 
plies to a general time dependent Hamiltonian H{t) with 
eigenfunctions \m{t)) and eigenvalues E m (t). The time 
derivative of \m(t)) should be calculated according to 



dt\m(t))=J2-^-\m(t)) x ^i, 



(22) 



with Ri the ith parameter. 

In the problem considered here, there are four param- 
eters: a, x, y, and z. Similar to the Ref. [5], we take the 
spatial mode functions to be identical as in Eq. (5). It 
is easy to see that the dark state Eq. (3) \D) as well as 
the other two eigenstates 



\B±) = 



1 



n 2 



(g\g ,i) 



i 



•%i,0»±-^|e,0>, (23) 



only depends on the parameter a and has nothing to do 
with r. So the numerator of Eq. (21) is independent of 
the atomic speed r. On the other hand, two of the three 
eigenvalues 



E±-- 



(24) 



do depend on r. Thus the denominator of Eq. (21) de- 
pends on the position of the atom r. Assuming that the 
time evolution of a(t) is uncorrelated with atomic motion 
r(t), the value of Eq. (21) may become large, especially 
in the regions where \E n (t) — E m {t)\ <C \(n(t)\d t \m(t)}\. 

Two regions require special attention: (1) the nodal 
planes perpendicular to the cavity axis due to the stand- 
ing wave term sin(fcz) = 0; and (2) the region away from 
the cavity axis due to the exponentially damped Gaus- 
sian term exp[— p 2 jui 2 (z)\. In this study, we focus on the 
level crossing dynamics, which happens mainly in the 
first region defined above. For the latest experiments, as 
in [1,2] where atoms are localized to a single well along 
the cavity axis, we may expect reduced nonadiabaticity 
because no actual level crossing occurs. 

Before attempting a comprehensive understanding of 
the problem, we will first discuss the qualitative picture of 
how adiabatic following is violated during the atomic mo- 
tion in this subsection. For simplicity, we will model the 
atomic motion as being simple one dimensional. From 
Eq. (13), the adiabatic condition Eq. (21) is just 



K(t) 
E±(t) 



< 1, 



(25) 



where \E±(t)\ — \x[ r (t)]\do\/^ + ko«(i)| 2 as defined be- 
fore. We see that E±(t) becomes sufficiently small (see 
Fig. 3) when the atom is near the nodal plane sin(kz) = 
0, where the condition Eq. (25) is easily violated. 




FIG. 3. The time dependent dressed state energies Eq/± 
with the time dependent pulse shape a(t) as in Eq. (14). 
We have used z(t) = vt with the speed of the atom being 
v — 0.14(A7) ~ 4 (m/s). a(t) is assumed to begin increasing 
at time t — when the atom is located at the peak of the 
cavity field where | sin(fcz)| = 1. 



After elementary substitutions, we find 



K(t) 



E±(t) 



fit) 



with 



/(*) = 



| sin(fcz)| 
a{t) 



v^oCi + K*)! 2 ) 3 / 2 

In Fig. 4, we have graphed the function f(t). 



(26) 



(27) 



0.04F 



0.03 



£-0.02 



0.01 




FIG. 4. The function f(t) of Eq. (27). 
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The two prominent features as in Fig. 4 at times 
of 3.3/7 an d 16. 7/7 correspond to the instants when 
the Gaussian shaped pulse gives rise to largest shape 
changes. Apparently, f(t) is significant at t = 3.3/7 an d 
t = 16. 7/7. From Eq. (26), we note if the atomic position 
is near a nodal plane when | sin(fcz)| is small, | e±(1) 1 
comes large and the adiabatic condition (25) can become 
severely violated at these instants. This adiabatic break- 
down is less serious at t = 16. 7/7 when the internal state 
transfer protocol is almost completed, but detrimental at 
t = 3.3/7, near the beginning of the process. 



III. LANDAU-ZENER TRANSITIONS 



( j_ =e - i J*e(t')dt'/n KCo 

= e - 2l l> (t ' )dt ' /h C + , (32) 
and C±(—u) = 0, we find 

c = r* e -2if*' e(t")dt"/h£ + ^ dt , 
J —u 

= e -2itf e(t")dt"/n c+ 

+2i f e- 2l Io £{t " )dt " /h e(t')C + (t')dt'/h. (33) 
Substituting Eq. (33) to (31), we find 



As was shown from the previous discussions, the appli- 
cation of the adiabatic approximation can potentially fail 
near a nodal plane where x(r) = 0. In this section, we 
hope to analytically investigate the transition from the 
dark state to bright states when the atom goes across a 
nodal plane. 



A. One-dimensional motion along the cavity axis 

We first deal with the one dimension case of atomic 
motion along the cavity axis. Assuming 



X (r) = sin(kz), 



(28) 



and take the atomic motion according to z = vt, our 
problem is to solve Eq. (12), with the initial conditions 



|0>(-«)l = 1. 
C±(-u) = 0. 



(29) 



More specifically, we in fact only wish to solve for |Co (u) \ 2 
and |C±(w)| 2 after one nodal crossing. 

As the lowest order approximation, we assume a(t) 
and K(t) to be constants in the domain kz £ 
[—u,u] and we also assume that the energy e{t)/h — 
sin(kz)go\/l + \a(t)\ 2 can be approximated as a linear 
function of time (see appendix B) in this domain 



e(t)/H w kvtg y/l + \a(t) | 2 , 



(30) 



where v is a constant atomic speed along the cavity axis. 
The second equation of Eq. (12) can be rewritten as 

= ie{t)C+/h + e l Io ^')dt'/n R( j o 

= ie(t)C+/h~e l f» £{t ' )dt ' /h K 

(e~ l Jo £{t ' )dt ' /n KC+ + e l> e(t ' )dt ' /h KC-) 
= ie(t)C+/h - K 2 C+ - e 2l $l e{t ' )dt ' /h K 2 C-. (31) 

Using 



—C + - ie(t)C/h - 2K 2 C+ 



+2iK 2 f e 2t l t ' £(t " )dt " /h e(t')C + (t')dt'/h. (34) 

J — u 

The last term is rapidly oscillating and within the lowest 
order approximation, it can be neglected (see appendix 
C). Then the equations of C+ and Co become 



dt 2 



C+=ie{t)C + /h-2K 2 C + , 



and 



C + =e l J« E(t ' )dt ' /h KC . 



(35) 



(36) 



Using the transformation 

d = (-i)-Lco, 



c + = e 



-i f*e(t')dt'/h 



C+, 



(37) 



we find that Eq. (35) is equivalent to a problem described 
by an effective Hamiltonian as below 



*4 ( do 

dt \ c+ 



( -y/2hK \(do\ ( ^ 
{-V2hK e ){c + )' (38) 



with the initial condition 



do(-u) = (-i) 
c+(-u) = 0. 



V2' 



(39) 



Nonadiabatic effect induced transitions mainly occur 
within the domain kz £ [—u,u], i.e. u is chosen such 
that beyond this domain, adiabatic condition Eq. (21) 
is well satisfied. In the end, as we will see later that 
our result is independent of the choice of u. Since e(t) is 
a linear function of t, this problem is exactly the same 
one as discussed in the original Zener's paper [10]. At 
the edge of the domain kz £ [—u,u], \K\ <C |e| and the 
initial conditions can be adiabatically maintained to the 
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domain (—00,00), just as Zener has done. Noting the 
normalization condition 



I do 



co 



and using Zener's solution [10], we find that 



MoMI 2 = 

which leads to 

|c + M| 2 

with 

h(t) 



■ exp 



-2tt- 



2K 2 



\kvg y/T 



(40) 



(41) 



|do(oo)| s 



1 — exp 



2t h(t) \ 
\kvgo\J 



2K 2 



d 2 (i) 



| a )3)B/3 



(42) 



(43) 



This constitutes the main result of our paper. We note 
the dark state probability after crossing a node becomes 



P 



1 - 
exp 



C + (oo)| 2 -|C_(oo)| 2 
h(t) 



-2tt 



\kvg Q \ 



(44) 



This leads to the conclusion that the larger the atomic 
speed v z is, the more reasonable it is to adopt the approx- 
imation of taking a{t) and K(t) as constants and neglect 
the term 2iK 2 f_ u exp[2i f* e{t")dt" /h]e{t')C+{t')dt' /%. 
Furthermore as we shall see in the next section, even in 
the limit of a small v z , the transition probability (44) 
as given by the analytic Landau-Zener method also com- 
pares well with results from numerical simulations. 



B. 3-dimensional motion of atoms 

In the limit as considered presently when atomic mo- 
tion is predetermined, a full 3-dimensional center of mass 
motion of the atom can be discussed without much fur- 
ther complications. Essentially, it is the component of 
the atomic velocity along the cavity axis direction that is 
involved in the level crossing dynamics, motion in the or- 
thogonal directions only causes the crossing to be at dif- 
ferent radial locations, thus different level spacing char- 
acteristics. 



Unless otherwise noted, the parameter go — (27r)50 
(MHz) is used, corresponds to the ~g at the end of Sect. 
II of Ref. [5]. The atomic mass for Cs is M « 133 x 
1.67 x 10~ 27 (kg). 

We assume the pump shape is given by Eq. (14) where 
a = 30, T = 3 x 10~ 7 (s) ~ IO7 (so that the total 
operation occurs within 2Xb = 2O/7, as typical for the 
optimal STIRAP process), and t\y — To/ 3. Therefore 
te(0,2T ). 

From the approximate Landau-Zener result Eq. (44), 
it is easy to see that the atom's dark state survival prob- 
ability P after crossing a node is closely related to the 
function h(t) of Eq. (43). The larger is h(t), the smaller 
is P. We show the time dependence of h(t) in Fig. 5. 
It resembles the function f(t) as shown before in Fig. 4. 
The duration of the state transfer protocol is taken to 
be 2O/7 so that the (unavoidable) maximums of h(t) are 
clearly displayed. 
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FIG. 5. The function h{t). 

It is easy to see that at times near i = 3.3/7 or < = 
16. 7/7, h(t) becomes rather large. If an atom crosses 
a node at these instants, the transition probability to 
other states may become significant. This corresponds 
to the qualitative picture of the adiabatic breakdown as 
mentioned in the last section. 

As a simple example, we assume the atom moves with 
a constant speed v. At time t = when the state transfer 
protocol begins, the atom is at the peak of the cavity field 
where | sin(fcz)| = 1. In Fig. 6 we present the results for 
the dark state survival probability P (as a function of v) 
at t = 2O/7, after the internal state transfer protocol has 
been completed. 



IV. RESULTS AND DISCUSSIONS 



In this section, we investigate the dark state survival 
probability by comparing the analytic result (44) under 
the Landau-Zener approximation with numerical solu- 
tions of Eq. (12). 
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FIG. 6. The dark sate survival probability P after crossing 
a nodal plane when the atomic motion corresponds to a con- 
stant speed along the cavity axis. The solid line comes from 
the numerical simulation of the nodal crossing dynamics by 
solving the Eqs. (12) for Cj(t), while the dashed line is the 
prediction of our approximate analytic result Eq. (44) from 
the Landau- Zener theory. 

The two prominent features of small P valleys can be 
easily understood. They correspond respectively to the 
crossing of a nodal plane at instants when h(t) is large as 
in Fig. 5, by slow and fast moving atoms. If the atomic 
speed is small enough, e.g. when v < 0.01 (A7) = 0.28 
(m/s), P remains essentially unity. This is because with 
this speed, the atom can't arrive at the nearest node 
before the state transfer protocol is completed. When v 
is increased to 0.017(A7) = 0.47 (m/s), the first small 
P valley shows up, corresponding to the atom arriving 
at the nodal plane at about t = 16. 7/7 when h(t) is 
significant (near its second peak). The second small P 
valley corresponds to the atomic speed of 0.049(A7) = 
1.36 (m/s), when the atom arrives at the nodal plane at 
about t = 3.3/7 when h(t) is around its first temporal 
peak. As we have analyzed in the previous section, this 
second peak corresponds to the beginning of the state 
transfer protocol. If the atom leaves the dark state at 
this time, the whole operation will be destroyed. 

We also note that the minimum of P near the sec- 
ond small P valley when v ~ 0.049(A7) = 1.36 (m/s) 
is larger than the minimum of the first small P val- 
ley of v ~ 0.017(A7) = 0.4 (m/s). This can be eas- 
ily explained according to the analytic result Eq. (44), 
P = exp[— 2nh(t)/\kvgo\] ) which shows that for the same 
value of h(t), P is larger for larger v. When the atomic 
speed is larger than v ~ 0.049(A7) = 1.36 (m/s), the 
atom will cross more than one node during the oper- 
ation time, and the dark state survival probability be- 
comes even smaller. We can approximate in this case 
P f» Y\i=i Pi where n is the number of nodes crossed by 
the atom and Pi is the probability for the atom to re- 



main in the dark state after crossing the ith node. This 
constitutes an excellent approximation when n is small. 

In practice, the atom may be trapped in an additional 
potential, e.g. takes a harmonic motion instead of a 
straight line. In the optimal scenario when the center 
of the harmonic trap overlaps the peak of the cavity field 
standing wave, and when the operation starts at the in- 
stant when the atom is located at the trap center, the 
corresponding results for this case is presented in Fig. 7, 
where we have further assumed a typical trap frequency 
ujt ~ 1-32 (MHz). We note that the atom's final dark 
state survival probability is related to its initial speed v 
as well. 



v (m/s) 

0J5 1 4 1J5 2.0 




0.4 



0.2 



0.01 0.03 0.05 0.07 
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FIG. 7. The same as in Fig. 6 but for a predetermined 
oscillation of the atom as confined in an external harmonic 
trap. 

This figure resembles that of Fig. 6, with the main dif- 
ference being the minimum of P in the two valleys being 
smaller here. This is because in an harmonic motion, its 
speed at the nodal plane is smaller than its initial speed 
v. 

It is remarkable that despite of the approximations 
used in deriving the analytic Landau-Zener transition 
rate Eq. (44), it gives rise to results that show an over- 
all agreement with the fully numerical simulations of Eq. 
(12). This demonstrates convincingly that at least in the 
parameter regime being considered by us, our result Eq. 
(44) captures the complete physics involved in this model 
problem. 

A. the velocity dependence 

To gain some understanding of the effects due to the 
unavoidable momentum distribution of the atom, we as- 
sume here a one dimensional distribution (for the speed 
of the atomic center of mass) 
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f(v) 



A 



exp 



(Aw) 2 



(45) 



centered at a central velocity w and with a distribution 
width Aw = 1.2 x 10" 3 (A7) [0.035 (m/s)], or about 10 
times Cs recoil velocity. The normalization constant is 
given by 



/■oo 

A(v ,Av) = / exp 
Jo 



(Av) 2 



dv. 



(46) 



In the following, we consider vq £ [0,5.3] x 10~ 2 (A7) 
i.e. [0,1.5] (m/s). We note that for each w , the above 
distribution (45) is essentially bounded from above by 
-Umax ~ «o + 2Aw. When v = 5.3 x 10~ 2 (A7) = 1.44 
(m/s), we find v max ~ 5.5 x 10~ 2 (A7) = 1.5 (m/s) and 
iw2T /(A/2) = 2.21, i.e. for v Q < 5.3 x 10~ 2 (A7), the 
atom will cross at most three nodes. 

The dark state survival probability as a function of vq 
can then be approximately computed according to 



where 

B(v y 



4X2T 



exp 



4X2T 



exp 



4x2T 



= B(v ) 
A(vo)' 

{Av) 2 
{v - w ) 2 ' 



4x2T 

+ / exp 

3A 



4x2T 



1x2 T ( | 



{Av) 2 
(v - v ) 2 



exp 



4x2T 



(Av) 2 
(v - wq) 2 " 1 



(Av) 2 



(47) 

dv 

Pi(v)dv 

P 1 (v)P 2 (v)dv 

Pi(v)P 2 (v)P 3 (v)dv, 



and Pi(v), P 2 (v), and P$(v) are respectively the dark 
state survival probability after crossing the first, the sec- 
ond, and the third node. In above discussion, we have 
assumed that the atomic initial position is zq = A/4, 
where the cavity field has its maximal value. The P(vq) 
for other values of zq ^ can be obtained similarly. 

In Figure 8, we have presented the numerically com- 
puted dark state survival probability P as a function of 
wo for several different atomic initial position Zq. 
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FIG. 8. P as a function of vo for the initial atomic locations 
along cavity axis zq = A/4, A/3, and 5A/12. 



B. 3-dimensional atomic motion 

To complete this study, we present selective results for 
the 3-dimensional atomic motion in this subsection. We 
selected two different situations where the atom is ini- 
tially at the anti-nodal point of the cavity field mode, 
and is taking a straight line motion that makes an angle 
of 30° or 60° with respect to the cavity axis. Not surpris- 
ingly, we again find excellent agreement with our analytic 
Landau-Zener result Eq. (44), applied appropriately as 
discussed earlier with the velocity component along the 
cavity axis being used to parameterize level crossing, es- 
sentially the same as the case of the 1 dimensional model 
considered earlier. 
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FIG. 9. The dark sate survival probability P after crossing 
a nodal plane when the atomic motion corresponds to a con- 
stant speed along the direction 30° off the cavity axis. The 
solid line comes from the numerical simulation of the nodal 
crossing dynamics by solving the Eqs. (12) for Cj(t), while 
the dashed line is the prediction of our approximate analytic 
result Eq. (44) from the Landau-Zener theory. 



To be sure of the adiabatic following of the dark state, 
one needs to assure at all times 



(27T) 



h(t) 
\kvg \ 



<1, 



(48) 



and the number of nodes crossed is small. 
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FIG. 10. The same as in Fig. 9 but for atomic motion 60° 
off the cavity axis. Only one small P valley shows up within 
the velocity range because of the large angle off the cavity 
axis. 



V. SUMMARY 

In conclusion, we have studied the nonadiabatic mo- 
tional effect of a three-level A-type atom Raman coupled 
to the standing wave quantum field of a high Q opti- 
cal cavity and an external pump field sharing the same 
spatial profile. 

First, making use of the Landau-Zener approximation 
to the crossing of a nodal plane by the atom, we have 
derived an analytic formula describing the survival prob- 
ability for the atom to stay in the so-called motional in- 
sensitive dark state. Surprisingly, our numerical results 
show that the approximation is remarkably good within 
current experimental parameters, thus can be used to 
guide the experimental implementation of the motional 
insensitive protocol [5]. 

Second, we find that the nonadiabatic motional effects 
is essentially connected with the dimensionless parameter 
v x (2O/7), the distance the atom (with center of mass 
velocity v) travels during the state transfer protocol of 
~ 2O/7. If this distance becomes a significant fraction 
of A, i.e. v x (20/7) > A/4, or v > A 7 /(48) = 0.577 
(m/s), then nonadiabatic effect will spoil the motional 
insensitive protocol in general, even if the atom is as- 
sumed to be located initially near the antinodal planes 
of sin(fcz) = ±1. 
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APPENDIX A: THE CASE OF A NONZERO 
DETUNING (A ^ 0) 

When A ^ 0, the Hamiltonian in the interaction pic- 
ture with the same basis |{|e, 0) , |s, 1) , \g, 0)} becomes 

/ A g fl \ 
H = [ g* . (Al) 

\n* J 

The three eigenvalues are 

E± = X - (A ± v/4 3 2 + 40 2 + A 2 ) , (A2) 

with the corresponding eigenstates \D) and |-E±). Clearly 
at nodal planes when sin(fcz) = 0, \E±\ take their mini- 
mal values 

|-E+|min = A, 

|£-|min = 0, (A3) 

as shown in Fig. 11, on inspecting of which leads to the 
following two comments. 

2i 7-^1 . 1 7-*r 1 




3. 1 
>* 

Oi 

5 0.5- 



CD 1 

Q -0.5 \ 



-1' ^ — ' ^ 1 

2 4 6 8 10 

kz 

FIG. 11. Similar to Fig. 2 but for A / 0. 
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First, irrespective of whether A = or A ^ 0, the 
dressed state energy levels cross at the nodal planes 
sin(fcz) = 0. 

Second, when A = 0, all three energy levels have the 
same value zero at the nodal planes, while for A ^ 0, 
only E- and E take zero values. There is a gap for 
E + whose width is A. Thus if A is large enough, the 
transition from dark state \D) to \E + ) can be avoided, 
but to state \EJ) remains because of the degeneracy at 
the crossing. The total transition probability again can 
be calculated theoretically using the previously adopted 
Landau-Zcncr approximation. 

To compute the transition probability for A ^ 0, we 
expand the state of the atom plus the field in terms of the 
eigenbasis Eqs. (3) and (8) of the system Hamiltonian for 
A = 0, which is now 



\D) 



1 



VT+Wt)\ 



-.(\ 9l ,0) - a(t)\g ,l)), 



\B ± ) = ^=(\B)±\e,0)), 



(A4) 



with the corresponding eigenvalues Eq. (7) reexpressed 
in this appendix as 



eo = 0, 

e±=±x[r(t)}\goWl + \a(t)\*. 
As before in Eq. (8), we have introduced 
1 



\B) 



+a(t)\g , 



(A5) 



(A6) 



and the parameter rn has been assumed to be unity. Then 
the system Hamiltonian including A is 

H = A|e,0)(e,0| + e+\B + )(B + \ + e.\B_)(B_\. (A7) 

Noting that 



|e,0) = -^(|B + )-|B_)), 



(A8) 



we rewrite Eq. (A7) as 



u= ie + + !) n.)(n.\ + L +|) /; /■• 



--{\B + ){B_\ + \B.){B + \ 



(A9) 



Expanding the quantum state as in Eq. (9) 

|*(t)) = C \D) + C+e- 1 ^ [e+{t ' )+A]dt ' \B+) 

+C-e- l $y-^+^ dt '\B_), (A10) 
we obtain the following equations 



C = 

C- = 
Since 



-(C+e- t Io [e{t ' )+A]dt ' + C_e" l /o he( *' )+A]dt V, 

C e l l> (t ' )+A]dt 'K 

+i | (C _ e <J>(*')^ + C+ ), 

C eJo l ~ e{t ' )+A]dt 'K 



+i-(C+e 



(All) 



(A12) 



C± = e ±i fo 2e{t ' )dt 'c^ 
we can integrate it to obtain 

C ± = j C T {t')e ±l 2 ^" )dt " dt> 

r~, ±i f* 2e(t')dt' 

~ C T e Jo , 

where we have , ne- 

glected the "small term" j C T (t')dexp[±i J*' 2e(t")dt"} 
(as in Sect. Ill and appendix C). Then according to Eq. 
(All), we have 



(A13) 



C + =C e i f> it ' )+A]dt 'K + l AC + , 
C.=C Q e i ^ t ' )+A]dt 'K + iAC. 



(A14) 



Assuming C± — e iAt £±, the above equation can be 
expressed as 

£t =Cve ±l 5l<^ dt ' K , (A15) 
which when coupled with the equation 

C = -(C + e- l f« Ht ' )+A]dt ' + C^I>^' )+A]dt ')K 
= -(£+e -< £ £(t ' ),tt '+e-e- i /o <^ dt ') K , (A16) 

is formally the same as equations for A = 0. Thus we 
obtain the same result 



ic | 2 = i-ie+i 2 -ic-i 2 
h(t) 



cxp 



-2tt 



\kvg \ 



(A17) 



APPENDIX B: THE LINEAR APPROXIMATION 

Within the discussion as in subsection II B, we mapped 
our level crossing problem into the well known problem 
of Landau-Zener transition. 

As was shown before, the adiabatic condition is 



fit) 
sin(fcz) 



«1, 



(Bl) 
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with the typical behavior for f(t) as shown in Fig. 4. 
Taking f(t) < 0.035, we see that within the domain of 
| sin(fcz)| < 0.7, we have 



/(*) 



sin(fcz) 



< 0.05, 



(B2) 



although still much less than 1. Thus we can define the 
domain | sin(£;,z)| < 0.7 as the domain of validity where 
the adiabatic condition is marginal. In this domain, the 
error of the linear approximation sin(/cz) ~ kz is about 
10%. 



APPENDIX C: THE SMALL TERM 

In this appendix, we provide the justification for the 
neglect of the second term of Eq. (34) . 

Given that initially the atom is in the dark state, we 
need C+(t') to be small in order to maintain adiabatic 
operation. Thus, we approximate 

c _ =e -2i f*'e(t")dt"/n c+ 

+2i [ e" 2l /o Eit " )dt " /n £ {t')C + {t')dt'/h, (CI) 



where u is the end of the time domain and assumed to 
satisfy | sin(A;z;u)| = 0.7. 

We first approximate the second term of Eq. (CI) 
according to 



2i f e - 2l fo e{t " )dt " /h e(t')C + (t')dt'/ 

J — U 

2/ / e' 2 ^o £{t " )dt " /h e(t')dt'C+(t)/h 



(C2) 



We note that 



S 



2i / e 



- 2i f° e{t " )dt " /h e{t')dt'/h 

-^t' 2 g ^l+\ a \ 2 kvt > go ^l + \ a \2 dt 



Ql 



e- tT T'dT 



(C3) 



where we have approximated a(t) as a constant and de- 
noted 



t' = nt'. 



(C4) 



S = 2 



\/ kvg y/l+\a\ 2 t 



(C5) 
(C6) 



where we have used sin -1 0.7 = 0.77. 

We take the worst case and use the value of a when a 
is maximum. This leads to 



a = 30 x exp 



0.55, 



fl x 10~ 7 -6 x 0.5 x 10- 7 ) 
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-14 



VT+H 5 = 1.14, 



(C7) 



and y kvgoy/l + \a\ 2 = 5.0 x lO 7 -^- If we now take 
t = s/(kv), we find 



yj kiig ^fl+\ a \ 2 Th 
-y/kvgoy/uWS^ 

e- lT '\'dT' 

5.69 



~ lT T'dT' 



ith 



x = y kvg Q y/l + \a\ 2 —. 



(C8) 
(C9) 

(C10) 



The oscillating behaviors of S for v = 3.5 x 10~ 2 (A 7 ) 
[1 (m/s)], 1.2 x 10- 2 (A 7 ) [0.35 (m/s)], and v = 3.5 x 
10~ 3 (A7) [0.1 (m/s)] are shown below in Figs. 12, 13, 
and 14. 
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FIG. 13. The same as in Fig. 12 but for v = 1.2x 10~ 2 (A7). 
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FIG. 14. The same as in Fig. 12 but for v = 3.5x 10~ 3 (A7). 
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